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We have used the very large JVAS/CLASS 8.4-GHz surveys of flat-spectrum radio 
sources to obtain a large, uniformly observed and calibrated, sample of radio source po- 
larizations. These are useful for many investigations of the properties of radio sources 
and the interstellar medium. We discuss comparisons with polarization measurements 
from this survey and from other large-scale surveys of polarization in flat-spectrum 
sources. 
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1 INTRODUCTION 

Polarization measurements are potentially important be- 
cause they provide information about magnetic fields and 
the media between source and observer. Radio emission from 
active extragalactic objects is generated by relativistic elec- 
trons spiralling around magnetic field lines, either close to 
the active centre or in the radio jets ejected from it. Radio 
waves passing through a magnetised plasma suffer Faraday 
rotation whose magnitude depends on the magnetic field and 
electron density along the path; this can occur within the 
synchrotron-producing region, in the interstellar medium of 
the emitting galaxy or in passage through our own Galaxy. It 
is also important to measure the polarization of large num- 
bers of discrete radio sources in order to remove their con- 
tribution from future measurements of the polarization of 
the cosmic microwave background radiation. 

Studies of radio polarization in extragalactic sources 
have yielded important results in the past, including sup- 
port for orientation models of active galaxies (Laing 1988; 
Garrington et al. 1988). They have also been used to investi- 
gate the physics of radio jets on kiloparsec scales, which can 
be polarized by tens of percent and whose polarization struc- 
ture can give information on magnetic fields and jet confine- 
ment (e.g. Laing et al. 2006; Leismann et al. 2005; Hardee 
2003; Gizani & Leahy 2003; Ferrari 1998; Laing 1996) down 
to parsec scales close to the central black hole (e.g. Lister 
& Homan 2005; Hughes 2005; Lyutikov, Pariev & Gabuzda 
2005; Fraix-Burnet 2002). 

The subject of this paper is the radio polarization of 
core-dominated, flat-spectrum, radio sources. At flux den- 



sity levels of >100 mjy at frequencies of a few GHz, these 
sources are predominantly quasars and BL Lac objects in 
which the bipolar jets are oriented close to the line of sight. 
Consequently, the radio emission from the approaching jets 
is enhanced by Doppler beaming, leading to associated ob- 
servable effects such as strong variability and apparent su- 
perluminal motion caused by the relativistic motion close to 
the line of sight. The integrated polarization in these flat- 
spectrum sources is typically a few per cent (e.g. Saikia, 
Kodali & Swamp 1985; Okudaira et al. 1993). The physics 
is somewhat complicated, since flat-spectrum "cores", when 
observed with milliarcsecond resolution, consist of multiple 
synchrotron components of different sizes and hence with 
different frequencies at which the emission becomes optically 
thick and below which the spectrum suffers a self-absorption 
turnover. The integrated rotation measure is usually quite 
low, often less than a few tens of rad m -2 (e.g. Rudnick & 
Jones 1983, but see also Zavala & Taylor 2004). 



Over and above these astrophysical inferences, occa- 
sional attempts have been made to investigate the hypothe- 
sis that observations of apparently ordered polarizations can 
be used to make cosmological inferences. Amongst these are 
chiral effects on the propagation of light (Nodland & Ral- 
ston 1997; see also Carroll & Field 1997), coupling of light 
pseudoscalars, such as axions, to the photons (Raffelt & 
Stodolsky 1988), and universal rotation (Birch 1982). More 
recently, Hutsemekers et al. (2005) have claimed alignments 
on large angular scales of the optical polarization position 
angles of quasars. Our data can be used to place constraints 
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on effects of this nature and these are discussed in more 
detail in Paper II of this series. 

In this paper we describe the data calibration and anal- 
ysis for the largest survey of polarization in compact ra- 
dio sources to date, namely the JVAS and CLASS surveys 
(Patnaik et al. 1992; Browne et al. 1998; Wilkinson et al. 
1998; Myers et al. 2003; Browne et al. 2003). We discuss 
the steps that need to be taken to ensure clean polarization 
measurements with low systematics in these large samples. 
In subsequent papers we address uses of this sample. In Pa- 
per II we discuss cosmological alignments prompted by the 
claims of Hutsemekers et al. (2005) of extreme-scale align- 
ments of quasar optical polarization vectors. In Paper III we 
will present further polarization observations of a subset of 
the JVAS/CLASS sources. 



2 OBSERVATIONS AND DATA ANALYSIS 

2.1 The JVAS/CLASS surveys 

The prime motivation for JVAS (the Jodrell-VLA Astromet- 
ric Survey) was to identify sources suitable for use as inter- 
ferometer phase calibrators and to measure accurate posi- 
tions for them. A secondary aim was to look for any sources 
that might have been gravitationally lensed. For both these 
objectives, pre-selecting sources which one might expect a 
priori to have compact radio structures was desirable. For 
this reason only those sources with known flat radio-spectra 
(S^.-Pghz > —0.5) were selected. Spectral selection for JVAS 
was initially carried out using the 5-GHz 87GB survey (Gre- 
gory & Condon 1991) and the 1.4-GHz Green Bank survey 
(White & Becker 1992). Each JVAS source was then ob- 
served for 2 min with the VLA in its A-configuration at a 
frequency of 8.4 GHz. A full description of the sample selec- 
tion process and the observational details can be found in 
Patnaik et al. (1992). In three separate observing sessions, 
each separated by around 15 months, 2720 sources stronger 
than ~200 mjy at 5 GHz were successfully observed in the 
region of sky with declination > 0° and |b|> 2.5°. The ob- 
servations in each session were typically divided into blocks 
of time lasting for >12 hr, each of which was calibrated and 
analysed separately. 

The CLASS survey (Cosmic Lens All-Sky Survey) con- 
sisted of 30-second snapshot observations at 8.4 GHz of all 
sources with a 5-GHz flux density of >30mJy which had not 
already been observed as part of JVAS. In addition, some 
sources below 30 mjy were observed, down to the GB6 sur- 
vey limit of ~20 mjy and using newer surveys (GB6, Gre- 
gory et al. 1996; WENSS, Rengelink et al. 1997; NVSS, Con- 
don et al. 1998). Full observing details and a more complete 
description of sample selection can be found in Myers et al. 
(2003). Four observing sessions in 1994, 1995, 1998 and 1999 
covered the whole sky between declinations of 0° and 70°, 
with the more southerly regions being observed later in the 
programme. 13783 sources were observed, making a total of 
16503 sources for the whole JVAS/CLASS programme. 

2.2 Calibration 

The original purpose of the CLASS survey was to identify 
gravitational lens systems. The analysis was done using a 



difmap script (Shepherd 1997) for automatic mapping and 
identification of candidates. For the present work we have re- 
edited, recalibrated and re-imaged the datasets as described 
below. This is because we have found that bad data on in- 
dividual telescopes or calibration errors which do not cause 
major problems for lens identification can become serious 
impediments to accurate polarization calibration and map- 
ping. All analysis was done using the NRAO Astronomical 
Image Processing System (aips). Because we want to be able 
to perform statistical analysis on large samples of sources 
with low levels of polarization, the polarization calibration 
is crucial, and we therefore describe it in some detail. Er- 
rors in the antenna polarization calibration at one epoch of 
observations will show up as correlated polarizations with 
similar position angles in the patch of sky observed at that 
epoch. 

After careful editing of bad data, phase and amplitude 
solutions were obtained using point-source phase calibrators 
which were typically observed every 15-20 minutes during 
the CLASS observations, and more frequently in the case 
of JVAS. These solutions were inspected, edited where nec- 
essary to remove discrepant points, and interpolated before 
application to the data on the target sources. 

2.2.1 Instrumental polarization terms 

Instrumental polarization terms were then estimated us- 
ing the AIPS task PCAL. For this process it is necessary to 
find a calibrator which either has a known polarization, or 
has an unknown polarization but which has been observed 
at a range of parallactic angles; this allows its polariza- 
tion at the time of observation to be calculated. For the 
majority of epochs of CLASS data the unpolarized source 
3C84 (J0319+415) was observed. Aller, Aller & Hughes 
(2003) present monitoring results for 3C84 over 17 years 
and quote a 14. 5-GHz polarization of 0.12%±0.01%, which 
is well within the error to which we are able to determine 
instrumental terms. We have therefore assumed 3C84 to be 
completely unpolarized. In some cases we have used OQ208 
as a zero-polarization calibrator which has been observed on 
numerous occasions with the WSRT (A.G. de Bruyn, pri- 
vate communication). In a few cases, particularly in many 
of the CLASS observations from 1994, we do not have a 
suitable instrumental term calibrator and have therefore ex- 
cluded these observations from our analysis. This results in 
a sparse sampling of some of the sky, in particular for areas 
north of declination 45° . 

For one observing session in 1991, we have observa- 
tions of both the unpolarized calibrator OQ208 and obser- 
vations of a strong polarized calibrator, B1611+343, taken 
at a range of parallactic angles. In Figure 1 we show the 
instrumental terms derived from each of these two observa- 
tions. They are entirely consistent with each other giving us 
confidence in the use of OQ208 as an unpolarized calibra- 
tor. Figure 2 shows the resulting polarizations derived for 
150 target sources using the two methods. Little difference 
is evident. If significantly polarized sources are considered 
(polarized flux density >lmjy), the rms deviation in posi- 
tion angles between the two calibrations is 6°. We have also 
experimented with the use of the solution interval for the 
determination of instrumental term solutions (between 15s 
and 60s, a standard value of 30s being eventually chosen), as 
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well as choice of reference antenna, and find that the result is 
robust to these choices. In practice, the mean instrumental 
term for each telescope was normally between 1% and 2%. 
Excursions outside this range could nearly always be traced 
back to bad data in one or more telescopes which upsets the 
overall solution. 



2.2.2 Polarization angle determination 

After the instrumental terms were determined, the correc- 
tions were applied to observations of a polarized source of 
known position angle in order to determine the relative 
phase of the left-hand and right-hand polarization channels. 
This is necessary in order to rotate the phases for all base- 
lines within each polarization channel so that the average 
phases on the left-hand and right-hand channels are consis- 
tent with each other and also with the known linear polariza- 
tion angle of the position-angle calibrator. We used 3C286 
(position angle 33°) and 3C48 (position angle —70°) for this 
purpose (Myers & Taylor 2006). 

We have checked that the phase rotations applied to 
left- and right-hand data in each epoch within a season do 
not differ significantly, although some drift is expected re- 
sulting from system changes and which the phase rotation 
is designed to remove. In most cases this was true to within 
20° (r.m.s. deviation between epochs was 19° in 1995 and 
14° in 1998), apart from a significant change which occurred 
between 1998 April 5 and 1998 May 21. 

When the instrumental polarization term correction 
and the L-R phase rotation are both determined, a polariza- 
tion angle can be derived for each baseline from observations 
of the polarized calibrator. Ideally, there should be a small 
scatter between baselines of less than 10° in these measure- 
ments, giving a very small error in the overall polarization 
angle averaged over ~300 baselines. In some epochs, all 
from the CLASS part of the survey, this desirable level of 
scatter is exceeded by factors of 3-4. However, in such epochs 
where two polarized calibrators (3C48 and 3C286) are avail- 
able, we can determine the L-R phase corrections from one 
source and use them to attempt to recover the polarization 
angle for the other source. In these cases, and in other cases 
where two polarized calibrators are available, the angle is 
recovered to within 10° of that expected, suggesting that 
there is no major systematic problem in the use of data 
from these epochs and in particular that derived polariza- 
tion angles have systematic calibration errors of < 10°. In 
most of the survey, including all the JVAS observations and 
the CLASS observations from 1998 March and April, the 
reliability is likely to be significantly better than this. 

Table 1 shows the JVAS/CLASS epochs with the cali- 
brators used for each, together with the important polariza- 
tion parameters which have been derived for each epoch. 

2.3 Imaging and the derivation of position angles 

The calibrated datasets were mapped using an automatic 
procedure implemented using an AIPS runfile. A shift was 
applied to the phases of the u — v data in order to centre 
the observations on the brightest point source found in the 
earlier automated analysis described in Myers et al. (2003). 
Sources were not included if they were not detected in the 



previous analysis, or if a source was located but judged to 
be extended on the basis of angular sizes (defined as being 
fitted by a Gaussian FWHM of 0'.'5 or greater). An initial 
map of 12'.'8xl2'.'8 was produced using natural weighting. If 
the source had previously been detected by the Myers et al. 
(2003) analysis with a flux density of 40mJy or greater, the 
initial map was used as an input for phase self-calibration 
using clean components brighter than the first negative; in 
cases where the flux density is greater than 80mJy one it- 
eration of amplitude self-calibration was also applied, using 
the results from the initial phase self-calibration as an input 
model. For the purposes of visual inspection, final images 
in Stokes /, Q and U were then made using the CLEAN al- 
gorithm (Hogbom 1974) using the task AIPS IMAGR, and 
combined into polarization level and position angle maps. 

Polarization parameters were extracted using fits to the 
processed (and in the case of strong sources, self-calibrated) 
u — v Stokes Q and U data using the task AIPS uvfit. In 
order for this algorithm to converge successfully, the posi- 
tion of the centre of polarized flux is required. This was 
calculated by making a total-intensity map and using the 
brightest point in this map as the Q or U centre, on the 
assumption that these coincide. 

We did not use the value derived from applying the 
CLEAN algorithm to the Q and U data, as this can cause 
the measured flux densities to be decreased by CLEAN bias 
(e.g. Condon et al. 1998). The model-fitting procedure is 
valid provided that there is no extended polarized structure 
such as a polarized radio jet, in which case it provides an 
averaged value. 

We checked for the overall uniformity of the distribution 
of the model-fitted polarization position angles and found, 
as expected, that the position angles derived from direct fits 
to the data are statistically uniform in the range from 0° 
to 180°. In contrast, angles derived from CLEANed Q and 
U data show peaks at 0°, 45°, 90° and 135°, corresponding 
to CLEAN bias forcing the Q and U fluxes towards zero. 
This effect of CLEAN bias on polarization determinations 
is unexpected and potentially important; it will be discussed 
further in a future paper (Battye et al., in preparation). 

In order to keep only significant polarization detections, 
we imposed a limit of 1 mjy in polarized flux in most sub- 
sequent analysis, corresponding to approximately 4a; below 
these levels the derived polarized flux is increasingly affected 
by positive bias (e.g. Simmons & Stewart 1985). However, 
the polarization parameters for all sources are available in 
the electronic version of this paper. An extract from the first 
page of the data table is shown in Table 2 and includes the 
position of the source, the epoch of the data and the derived 
flux densities in Stokes I, Q and U together with the polar- 
ization percentage and position angle. The total number of 
sources in the catalogue with polarization measurements is 
12555 of the total JVAS/CLASS sample of 16503. Of these, 
4294 have significant detections of polarization, with a po- 
larized flux density >1 mjy. 
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Figure 1. Instrumental polarization terms per telescope for the observations taken in 1991 as part of the JVAS survey. The scales on 
the axes run from —3% to 3% of instrumental polarization. The length of each line gives the amplitude of the instrumental term; that 
is, the amplitude of polarization that this telescope would record for an unpolarized source. The angle gives the phase in such a way 
that a complete rotation on the diagram would rotate the polarization vector by 180°. The number attached to each arrow refers to the 
specific VLA antenna. On the left is the calibration achieved by using the polarized source 1611+343 observed at different parallactic 
angles, and on the right is the calibration made using the assumption that the source OQ208 is unpolarized. Note the similarity in the 
derived instrumental terms. 
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Figure 2. Comparison of derived polarization of 150 sources obtained using two different calibrations of instrumental polarization in 
observations from 1991 June 16. Axes represent Q and U as a percentage of Stokes /. Unfilled and filled circles represent calibrations 
using the polarized source 1611+343 and the unpolarized source OQ208. Larger symbols represent significant detections of polarization 
(>1 mjy in polarized flux density). The right-hand panel is a blowup of the inner region of the left-hand panel, indicated by a dotted 
line in the left-hand panel. 



2.4 External comparisons with other data. 

Do the CLASS* polarization position angles agree with 
other measurements? We have made comparisons, some di- 



* In this and subsequent sections we refer to the combined 
JVAS/CLASS surveys simply as "CLASS" . 



rect and some indirect, that convince us that our CLASS 
polarization angles are reliable. Okudaira et al. (1993) have 
used the Nobeyama 45-m telescope to measure the polar- 
izations of 99 flat spectrum radio sources at 8.4 GHz and 
their sample has 76 objects in common with ours. We plot 
the histogram of position-angle differences in Fig. 3. There 
is a clear peak around zero degrees, with a rms dispersion of 
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Calibrator from 950814 used 



3C286 agrees to <10% in R-L off- 
set despite large scatter 



Calibrator from 980314A used 



Calibrator from 980521 used 



Table 1. Observing epochs for JVAS/CLASS. For each epoch we list the approximate sky area principally covered by the observation. 
The zero-polarization calibrator observed to calibrate the instrumental polarization terms is listed together with the derived instrumental 
polarization in columns 4 and 5; in most cases these can be brought below 1.5% with careful editing of bad data. The amplitudes of 
these terms arc also given by the lengths of the vectors in Fig. 1. Column 6 lists the polarized calibrator used to determine L-R phase 
offsets as described in the text. The scatter between measurements of polarization angle for different baselines, derived from applying 
the instrumental polarizations to the polarized calibrator, varies widely between epochs. 



±40° . Some spread is expected because many of the sources 
are known to be extremely variable both in their total in- 
tensity and polarization properties, position angles (e.g. BL 
Lac, 3C345, 3C454.3, etc.). In addition, the angular reso- 
lutions of the two sets of observations are very different. 
Resolution may make a difference because emission from 
optically thin, kiloparsec-scale jets, is generally more highly 
polarized, and sometimes at a different position angle, com- 
pared to that from the compact cores (e.g. Saikia & Salter 
1988). 

Another test is to compare the CLASS position angles 
with the intrinsic position angles, corrected for Faraday ro- 
tation, determined by Broten et al. (1988). Typical rotation 
measures arising from a combination of Galactic and inter- 
nal Faraday rotation are ~20 rad m~ 2 which means that 
most of the CLASS measured angles should be close to the 
intrinsic ones. The difference between the CLASS angles and 
those from Broten et al. (1998) for the 67 objects in com- 
mon are also shown in Fig. 3. The spread is about the same 
(~ 40°) as that of the Okudaira et al. (1993) sample but 
there is a peak which a formal average of the data shows to 
be centred around ( — 7 ± 5)°. 

A typical rotation measure of ~ 20 rad m~ 2 will pro- 
duce a rotation of ~ 50° at a wavelength of 21 cm which, if 
left uncorrected, would effectively smear out any record of 
the intrinsic position angle. There are, however, areas of sky 
at high Galactic latitudes where average rotation measures 



are significantly less than 20 rad m -2 and hence one would 
expect to be able to see a correlation between the CLASS 
position angles and those measured at 21 cm in these re- 
gions. The NVSS survey, which was made at a wavelength 
of 21 cm, lists polarization information on all the sources 
(Condon et al. 1998). In Fig. 3 we show the histogram of 
position angle differences between CLASS and NVSS for the 
1424 common sources in the region |b|> 30° and which have 
NVSS polarized flux density >2 mjy and CLASS polarized 
flux density > 1 mjy. As expected there is a significant peak 
near zero degrees, although the mean value of the data is 
actually located at (3.9±1.2)°. 

The offset in the peak from 0° in the CLASS-NVSS po- 
larization angle difference is almost certainly due in part 
to large-scale rotation measure structure in the Galaxy. 
Simard-Normandin & Kronberg (1980) present RM mea- 
surements over the whole sky, and identify regions of 60°-90° 
in extent with correlated Galactic RM. To demonstrate that 
the observed offset in the peak is probably due to Galactic 
Faraday rotation, we have selected sources in regions found 
to have significant positive RM (defined by 90° < I < 165°, 
b > 0°; Simard-Normandin & Kronberg 1980) and signifi- 
cant negative RM (defined by 5 > 0°, I < 130° and b < 0°, 
corresponding to Simard-Normandin & Kronberg's region 
A). Plotting these separately (Fig. 4) clearly shows that the 
polarization differences are offset in different directions in 
these two regions of sky. 
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Table 2. The first page of the polarization angle measurements with the automated analysis as described in the text. The first six 
columns are the position in RA (3 columns) and declination (3 columns) of the major component as derived from the CLASS analysis 
(Myers ct al. 2003). Remaining columns contain measurements derived from the direct fits to the u — v data. Flux densities in millijanskys 
follow for total flux (Stokes /), error in /, Stokes Q, error in Q, Stokes U and error in U. The final column is the derived position angle 
in degrees (i tan -1 (U/Q)). 



Finally we emphasize that seeing any peak in the 
CLASS-NVSS position angle histogram implies that the av- 
erage integrated rotation measures of the sources in the sam- 
ple must be < 20 rad m~ 2 . This is somewhat surprising 
since in these flat-spectrum sources it is believed that most 
of the emission originates on parsec-scale structures, and it 
is now well established that the core regions detected in cm- 
wavelength VLBI observations typically have rotation mea- 
sures of several hundreds of radm~ (e.g. Zavala & Taylor 
2004). Indeed, the core rotation measures revealed by 7mm- 
2cm VLBI observations sometimes approach 1000 radm~ 2 
or more (Mutel et al. 2005; Gabuzda et al. 2006). 

We have also performed an indirect test. It is known 
both statistically, and from detailed mapping of individ- 
ual radio sources, that there is a correlation between the 
structure of radio jets and the magnetic field direction in- 
ferred from radio polarization measurements. In particular 
high-luminosity quasar jets tend to have dominant longitu- 
dinal magnetic fields (Bridle & Perley 1984; Cawthorne et 
al. 1993). Therefore we might expect to see a correlation be- 
tween jet position angles and polarization angles. We have 
compiled a list of 157 sources which have VLBI elongation 
directions given in the tables of Taylor et al. (1994) and 
Henstock et al. (1995). In Figure 5 we show the histogram 
of the difference between structural and polarization posi- 
tion angles for these sources in common between the two 
samples. There is a clear excess around 90 degrees differ- 
ence, as expected if magnetic field direction aligns with that 
of the radio jet. 



3 CONCLUSIONS 

We have presented polarization determinations for the 76% 
of CLASS sources for which reliable polarization calibration 
is possible. Although systematic effects are important in de- 
termining low levels of polarization, we have quantified the 
systematic errors involved and compared the data with other 
surveys. Agreement is obtained which convinces us that the 
polarization angles are reliable to within a few degrees for 
sources with a significant (>1 mjy) polarized flux density. 
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Figure 4. CLASS-NVSS polarization position angle differences for two regions of sky. On the left is the region 90° < / < 165° , b > 0° , 
found to have significant positive RM by Simard-Normandin & Kronberg (1980). On the right is the region of significant negative RM, 
<5 > 0°, I < 130° and b < 0°. 
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Figure 5. Histogram of position angle offsets, in degrees, be- 
tween CLASS polarization vectors and VLBI structure from the 
observations of Henstock et al. (1995) and Taylor et al. (1994). 
Note the peak at about 90°. 
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Figure 3. Histogram of position angle differences, in degrees, be- 
tween CLASS and other surveys. From top to bottom: Okudaira 
et al. (1993), Broten et al. (1988) and NVSS (Condon et al. 1998). 
In the case of NVSS we have required sources with significant po- 
larization in each survey (>1 mjy for CLASS and >2 mjy for 
NVSS). 
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